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Abstract. - The two-component Fermi gas in a double- well trap is studied using the density func- 
tional theory and the density profile of each component is calculated within the Thomas-Fermi 
approximation. We show that the two components are spatially separate in the two wells once the 
repulsive interaction exceeds the Stoner point, signaling the occurrence of the ferromagnetic tran- 
sition. Therefore, the double-well trap helps to explore itinerant ferromagnetism in atomic Fermi 
gases, since the spontaneous separation can be examined by measuring component populations in 
one well. 
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Introduction. -Itinerant ferromagnetism has long 
been one of the central research topics in condensed matter 
physics [I] . Recently, developments of laser trapping and 
coohng of atomic gases with multiple components have 
greatly renewed and broadened the interest in this area, 
and the new research interest is mainly expressed in two 
aspects. 

First, the *^Rb atoms provide an opportunity to study 
the itinerant-boson ferromagnetism [21 ■ It is the 
itinerant-electron ferromagnetism that has been inten- 
sively studied in the context of condensed matter physics. 
That is a typical Fermi system. The ®''Rb gas comes as the 
first example of the ferromagnetic Bose system. Motivated 
by the experimental achievement, the ground state prop- 
erties [His] and thermodynamics [HIT] of the Bose ferro- 
magnet have been studied theoretically. It was shown that 
the Bose gas is much easier to exhibit ferromagnetism than 
the Fermi gas: In the former case, the ferromagnetic (FM) 
transition temperature is never below the Bose-Einstein 
condensation temperature regardless of the magnitude of 
the ferromagnetic coupling [6], whereas in the latter case 
the ferromagnetism can not be present unless the ferro- 
magnetic coupling exceeds the Stoner point [T]. 

Second, cold atomic Fermi gas can be used to simu- 
late mechanism of the itinerant-fermion ferromagnetism, 
a long unsolved question in condensed matter physics. So 
far, the Stoner model gives us a qualitative description of 
itinerant ferromagnetism [1 . According to this model, the 
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electron system can lower its total energy by spin polariza- 
tion when the decrease of interacting energy is larger than 
the increase of kinetic energy due to Pauli principle. Since 
the interaction between atoms is tunable by Feshbach res- 
onance [5], people can examine whether such transition 
occurs when the repulsive interaction becomes stronger. 

Recently an experimental group from MIT claimed their 
realization of ferromagnetic phase of Fermi gases in an 
equally populated mixture of ^Li atoms in the lowest 
two hyperfine states [5]. They addressed the observa- 
tion of non-monotonic behavior for increasing repulsive 
interactions, which implies the occurrence of itinerant FM 
transition through comparing with the Stoner's ferromag- 
netic mean-field theory. Previously, Duine and MacDon- 
ald had already investigated features of the FM transition 
based on second-order perturbation theory [TU]. LeBlanc 
et al. studied observable experimental signatures of the 
FM transition theoretically within a local density approx- 
imation [TT]. More elaborate theoretical studies beyond 
the mean-field approximation showed that the FM transi- 
tion could take place at a weaker interaction strength [T^] . 
Moreover, various related problems, such as textured mag- 
netization [13], pairing instability [14], spin fluctuations 
[15j , and population imbalance |16j , were investigated. On 
the other hand, the existence of a FM transition was ques- 
tioned by Zhai, who suggested in a phenomenological way 
that a correlated state without ferromagnetism could also 
cause the similar features in the experiment [17 . 

Direct evidence of the FM transition in atomic gases 
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might be the formation of magnetic domains, or the phase 
separation of different spin components. Amoruso et al. 
[18] and Salasnich et al. [19] calculated the density pro- 
file of two-component Fermi gases with conserved particle 
numbers and demonstrated the spatial symmetry break- 
ing of the system. This means the formation of a sort of 
domain-structures. Furthermore, Sogo et al. indicated the 
similar phenomenon in an open system with variable par- 
ticle numbers [20] . The MIT experiment failed in viewing 
the domain structure. 

It is not an easy job to observe domain-structures in 
atomic gases directly. As Jo et al. claimed, the lifetime 
of domains might be too short and their size might be 
too small to be observed [21] . Phase-contrast imaging has 
been successfully used to measure magnetization in situ 
in ®''Rb condensates [3], but in situ detection for atomic 
Fermi gases has not been reported as yet. Some recent 
theoretical proposals seek to verify the FM transition in- 
directly, e.g., by detecting topological signatures [l^ and 
by calculating the spin drag relaxation rate [23] . 

Here we propose an alternative method to detect the 
spatial symmetry breaking in the ferromagnetic Fermi 
gas. This proposal involves trapping the two-component 
atomic gas in the double-well (DW) potential. We show 
that the two components are spatially separated and each 
component prefers staying in one trap if the FM transi- 
tion is induced by the strong repulsive interaction. This 
makes it convenient to determine the FM transition by 
measuring atoms in one well, in comparison to measuring 
magnetization in situ. 

This paper is organized as follows. In section II, we cal- 
culate the ground state energy of two-component Fermi 
gases in a DW trap using the density functional theory. 
The kinetic term we are using is due to the so called 
Thomas- Fermi approximation as used by Sogo et al. [20] . 
In section III, we show numerical results for the case with 
symmetric potential and equal particle numbers by ana- 
lyzing the Euler equation. In Section IV, we consider the 
asymmetric case by changing the potential depth and the 
particle numbers, respectively. The last section shows our 
conclusions. 

Basic model. -We now consider the ground state of 
two-component Fermi gases trapped in double wells at zero 
temperature. According to the density functional theory, 
the ground state energy of the system can be written as a 
functional of the densities: 
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where the first term on the right hand corresponds to the 
kinetic energy of the system known as Thomas-Fermi ap- 
proximation, and it has the similar form as the energy of 
uniform Fermi gases. The Thomas-Fermi approximation 
is considered a well approach near the center of the poten- 



tial trap for large total particle numbers [21]. The second 
term U{r)n„ is the potential energy of the system within 
local density approximation, and 
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describes the symmetric DW potential. Here x and p are 
the axial and radial coordinate. A and B are the constants 
related to the shape of the potential well, and cr = 1,2 
represents the two different component. gnin2 is the in- 
teracting energy, where the strength of coupling is given 
hy g = 4:TTh'^a/m with a the s-wave scattering length. The 
interaction between the same component is absent as a 
result of the Pauli exclusion principle. 

For simplicity, we transform the integral of equation (jlj 
into cylindrical coordinate and then convert the parame- 
ters to dimensionless forms, 
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Here f = ^/KJ muj is the oscillator length of the trap and 7 
is a variable relevant to the density profile and the particle 
number, which is determined by 
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where N„ is the number of each component and N„ is the 
reduced particle number. 

Then equation ([1]) can be simplified as 
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where U{x, p) = x^ + ce^^'-^' + p^ is the reduced potential 
with c, h are constants transformed from parameters A 
and B, respectively. 
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The ground state energy of the system can be obtained 
by minimize the energy functional. To ensure the con- 
servation of particle number of each component, it's nec- 
essary to introduce two Lagrange multiplier pi and p2, 
which just represent the reduced chemical potentials. The 
real chemical potential of the system /lo- reads 
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Fig. 1: Density profile along the axial coordinate x with radial 
coordinate p = 0. The reduced interaction strength a from (a) 
to (d) is 0, 0.9, 1, 1.05. The dotted-black curve illustrates the 
potential well. The solid-red and dashed-blue curves represent 
different components. 



The ground state energy has to fulfill the variational 
condition, 6{e — fiiNi — IJ.2N2) / Sn^ = 0. So we obtain the 
following Euler equations, 
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where Mcr=i.2{x, p) — fi^ ~ [x^ + ce 



p^). By solving 



the Euler equations and keeping the particle number of 
each component conserved, we can get the density profile 
of each component in the ground state. 

The symmetric case: with equal well depth and 
population. - In the case that the two-component 
fermions with equal population are confined in a sym- 
metric DW trap, it is reasonable to assume the chemical 
potentials of the two components equal, Hi = 1^2 = 
For an ideal Fermi system the densities of the two compo- 
nents will be the same everywhere in the DW trap. In the 
presence of weak interactions, the density will not change 
much. This corresponds to the symmetric solution of Eqs. 
pT|) and In addition, Eqs. (dJ) and ^ have the 

asymmetric solution, which means that density profiles of 
the two components become different. If the asymmet- 
ric solution is stable, the symmetry of the ground state is 
spontaneously broken. 

To minimize the ground state energy functional, one 
must make sure that d^e/dri^ > 0. This demands that 
the following inequality holds everywhere in the DW trap. 
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By analyzing Eqs. ((TT|) ^ P^ . the asymmetric solution is 
stable in case that the following condition is satisfied. 
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Fig. 2: Illustration of the density profile ni along the x and y 
coordinate with the interaction parameter a — 1.05. Numbers 
in the figure show the values of the density. 



There is only the symmetric solution when < M{x, p) < 
2!^, that is, ni = n2- Note that there is no stable so- 
lutions when M{x,p) > ^ according to Eq. ([13]), which 
means that the two components are completely separated, 
so we have n2 = and the equation for ni becomes: 
ni = M{x, p)2 . 

In the whole article we consider a gas composed of ^Li 
atoms with the trap frequency uj « (170 x 27r)IIz. First, 
the total number of two components is supposed to be A'^ = 
10^ (N^ = 0.5 X 10*5) and the DW potential is U{x,p) = 

-|-0.5e~^°^ + p^. Here the reduced particle number No- 
is normalized to 0.1 for convenience and 7 w 2.57x lO^m"^ 
correspondingly. The potential barrier is defined as Uq = 
\m,u?A, then C/q 15.35/i x KHz according to Eqs. 

Figure [T] portrays density profiles of the two compo- 
nents. The spatial width of the confined Fermi gas and 
the distance between the two wells are about 2 and 0.8. 
According to Eqs. (O, the corresponding physical values 
are about 119^m and 48^m, respectively. It indicates that 
the symmetry-breaking solution becomes stable as long as 
the interaction strength is large enough. The two com- 
ponents tend to spatially separated and each component 
prefers to gather in one trap. That is the ferromagnetic 
state of two-component Fermi gas in double wells. This 
state begins to appear at about the interaction strength 
ac ~ 0.9, corresponding to the critical scattering length 
flc = ic/7 ~ 6600ao, where ao = 0.529A is the Bohr ra- 
dius. As the interaction strength increases, the two com- 
ponents become more separated. Note that the compo- 
nents are separated in the inside region of each trap, but 
are still equally mixed in the outside region. Figure [2] il- 
lustrates the density profile of one component along the x 
and y coordinate with a = 1.05. 



p-3 



Junjun Xu et al. 



7000 




30 

U^(hxKHz) 

Fig. 3: The critical scattering length Uc/ao as a function of 
the particle number Na (a) and the barrier Uo (b). Here ao is 
the Bohr radius and h is the Planck constant. AfCT=i,2 ~ 0-5N , 
where A*' is the total number. 
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Fig. 4: Density profile along the axial coordinate x with radial 
coordinate p = with unequal particle numbers A'^i = 2N2 — 
0.2. The reduced interaction strength a from (a) to (d) is 0/0.5, 
0.8, 0.9, 0.94. The dotted-black curve illustrates the potential 
well. The solid-red and dashed-blue curves represent different 
components. 



LeBlanc et al. estimated that the critical value for 10^ 
^Li atoms confined in the single- well with the trap fre- 
quency uj « (170 X 27r)Hz is about 6400ao [H]. Therefore, 
trapping atoms in double wells does not cause much dif- 
ficulty in generating the FM transition, in comparison to 
in the single-well trap. On the contrary, it helps the two 
components to be separate in two wells, and thus produces 
benefits for the detection. 

Then we discuss the influence of the total particle num- 
ber and the potential barrier to the Stoner point. Fig- 
ure [3^ shows that the critical scattering length Oc de- 
creases apparently with the particle number increasing. 
Increasing the barrier Uq but keeping the particle number 
N„ = 0.5 X 10^ unchanged, the Stoner point is also low- 
ered, but not that significantly, as shown in Fig. It is 

worth noting that there is no significant change in kpUc 

1 

although flc varies, kp = [67r^no-(0)] ^ is the Fermi wave 
vector of the interacting gas with no-(O) the density at the 
center of each well. Our calculation shows kpUc ~ 1.57, 
which is consistent with the mean-field value in the uni- 
form case, kpCic — f ■ 

The asymmetric case. -To proceed, we consider the 
asymmetric cases caused either by the population imbal- 
ance, or by the difference of well depth. The density profile 
can still be obtained by numerically solving Eqs. (fTTj) and 
(fT2|). with the particle number of each component being 
conserved. The problem of fermion mixtures in a single 
well with population imbalance has already been investi- 
gated [16]. 

Figure |4] shows the results for the unequal-population 
gas in the symmetric DW trap. The ratio of total numbers 
A^i : A^2 is 2 : 1. Similar to the equal population case, 
the two components distribute symmetrically in the two 
wells when the interaction strength is relatively weak. As 



the interaction grows stronger, this symmetric distribution 
becomes unstable. The two components tend to repel each 
other and each component begins to occupy the center 
of one potential well, signaling the occurrence of the FM 
transition. 

Now look at a two-component Fermi gas with equal pop- 
ulation iVi = A'2 = 0.1 in an asymmetric DW trap. The 
trapping potential is chosen as U {x, p) = + 0.5e^^"^ + 
0.5sm(0.5a;) -I- . The third term in the right hand makes 
the trap asymmetric in the x-direction. As indicated in 
Fig. El the density of each component equals anywhere 
in the double trap in the weak interaction case, although 
they distribute asymmetric. As the interaction strength 
grows, the spatial separation of the two components oc- 
curs first in the deeper well, then in both wells. This 
phenomenon is easy to understand since the deeper well 
traps more atoms with larger density, so the interacting 
energy density is higher than that in the shallower well. 

At last, we illustrate the reduced chemical potential p, 
for both the symmetric and asymmetric cases in Fig. [6l It 
can be seen that the chemical potential monotonously in- 
creases with the interaction strength for all cases. For the 
symmetric-potential and equal-population case, the two 
components have the same chemical potential, even if they 
are spatially separated in the ferromagnetic state. In con- 
trast, for the asymmetric-potential and equal-population 
case, chemical potentials of the two components have the 
same value at first, and then become different when the 
two components start to separate (see the black arrow in 
Fig. E]). In the case of unequal-population, each compo- 
nent has different chemical potential. Comparing Fig. |4] 
with Fig. [5l it seems that the density profiles for the two 
asymmetric cases display some similar features when the 
interaction is strong enough, as shown in Figs.HJi and[5ji. 
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Fig. 5: Density profile along the axial coordinate x with ra- 
dial coordinate p = in the asymmetric double-well potential 
U{x,p) = + 0.56-^°="' + 0.5sm(0.5x) + . The reduced 
interaction strength a from (a) to (d) is 0, 0.8, 0.9, 1. The 
dotted-black curve illustrates the potential well. The solid-red 
and dashed-blue curves represent different components. 



Conclusion. -In conclusion, we have studied density 
profiles of two-component Fermi gases in double wells us- 
ing density functional theory. Both symmetric and asym- 
metric DW potentials have been taken into account and 
the unequal-population Fermi gas in the symmetric DW 
trap has also been discussed. For all the cases, the ob- 
tained results indicate that the two components can be 
spatially separated if the repulsive interaction becomes 
strong enough. This implies that the ferromagnet aligns in 
the z (population-imbalance) direction. In the ferromag- 
netic state, each component tends to gather in one well 
and thus its population dominates in this well. Therefore, 
the occurrence of the FM transition can be examined by 
measuring atomic populations in one well. It is much eas- 
ier than detecting magnetic domains in situ. 

Finally, we note that our calculations may not be di- 
rectly applicable when the ferromagnetic alignment is in 
an in-plane direction. Moreover, the present study can be 
extended in several directions. The critical point needs 
to be calculated more accurately. Calculations beyond 
mean-field theory and numerical simulations reveal a lower 
kpa « 0.8 — 1.1 in the single well [12]. It is natural to 
expect that it is also the case in the double well. Recent 
works also suggested that some nontrivial effects may play 
important role in understanding the FM transition, such 
as spatially modulated magnetism [T3] , pairing effect [2] , 
and spin fluctuations |15) . 
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Fig. 6: The reduced chemical potential /i versus the interaction 
parameter a for different cases. Here s is short for symmetric 
case, while ud and up correspond to unequal depth of potential 
well and unequal population. Note that fj, = pi = fj,2 for s case 
and the black arrow shows the transition point in the ud case. 
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